On the Long Distance Contribution to the B s — ► 77 Decay in the 

Effective Lagrangian Approach 

Wei Liu, Bin Zhang and Hanqing Zheng 

Department of Physics, Peking University, Beijing 100871, People's Republic of China 



'q^ ' Abstract 
On 
^ . 

We re-estimate the decay branching ratio of B s — > 27 through DfD~ and 
D* + D*~ intermediate states, in the effective Lagrangian approach. We find 
that the branching ratio does not exceed a few times 10 -7 , contrary to the 
result recently claimed in the literature. 
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The B s — > 77 decay is an intereting process which has been investigated by several 
groups fl|-|3[ within the context of pQCD. The branching ratio was found to be BTZ(B S — > 
77) ~ 3.8 x 10~ 7 for m t ~ 175GeV, to be compared with the present experimental upper 
limit BTZ(B S — > 77) < 1.48 x 1CT 4 (§]. It was also pointed out that the branching ratio 
can be substantially enhanced in some extensions of the standard model such as a generic 
2-Higgs doublet model ||. Provide that the short distance contribution is not large |||5],[7||, 
one may then hope the observation of such channel in future experiments as a signal of 
physics beyond standard model. However, as have been emphasized by Choudhury and 
Ellis [[|, a careful analysis to the long distance effects has to be made before one could use 
this channel as a probe to the physics beyond the standard model. Choudhury and Ellis 
have considered the contributions due to intermediate D s and D* states via the diagrams 
including loops of D s mesons alone, loops of D* mesons alone (see fig. [I]), and diagrams 
involving radiative D* — > D s + 7 transitions. Based upon an estimation on the absorptive 
part of the amplitude they find that the contribution due to D s alone is, 

Br(B s -> 2D S -> 77) ~ 2.9 x KT 8 , (1) 

which is insignificant as comparing with the pQCD short distance contribution. Diagrams 
involving the D* — > D s ^y transition are again found to be small. However they find that the 
contribution due to D* is much larger, 

Br(B s -> 2D* S -»■ 77) ~ 6.5 x , (2) 

which, if correct, is of great phenomenological interets as it will be easily seen in future 
hadron colliders. Because of its importance it is worthwhile to give a careful re-analysis to 
such a process. 
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FIG. 1. One-loop contribution to the B s — > 77 amplitude due to D s (D* s ) alone. 

We in the following use the same method as was adopted in ref. ||, i.e., a diagramatic 
analysis of fig. [1] based on an effective Lagrangian approach. Beside those having been done 
in ref. || we also estimate the real part contribution of the rescattering amplitude using 
dispersion relations. The Feynman rules we use for the calculation are derived from the 
dimension 4, U em (l) gauge invariant effective Lagrangian which are bilinear in D s fields: 



eff 



y»D- v - V u D-» ) + iefoD+F^D- , (3) 
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where <pv is an unknown parameter not constrained by gauge invariance. We have verified 
that the choice made in ref. || corresponds to <py = 0. In the case of elementary W + W~ r y 
interaction (fiv = 1 and in the case of effective p + p~7 or K* + K*~~f interactions the parameter 
4>v can be estimated from the extended Nambu-Jona-Lasinio model and is found to be 
slightly less than 1 || (Our definition of <pv is two times larger than Ref. ||). Since we 
did not find convincing reason to support taking (py = 0, we treat 0y as a free parameter 
ranging from to 1. 

The absorptive part contribution of fig. [I] can be written as, 



ImA = £A AlA2 /?T Al 



(4) 



where /3 = J 1 — 4m 2 /s is the kinematic factor, the sum is over the helicity indices Ai, A2 
of the intermediate states (for D* + and D*~ scattering) and A3, A 4 are the helicity indices 
of the outgoing photons; (Ai,A 2 ) can take (+,+), (— , — ) and (0,0), and (A 3 , A 4 ) can take 
(+, +) and (— , — ). In above Aq is the decay vertex obtained from factorization method, 

A = {Df(k)D;(q)\H wk \B s (p)) = -igf [(p 2 - k 2 )f+(q 2 ) + q 2 f-{q 2 )] , (5) 

for B s — > D+D~ , and, 

A = (D*+(k) D*-(q)\H wk \B 8 (p)) 

= 9f*eT£- [WV^Pah + i{A x {q 2 )v 2 r - 2A 2 (q 2 )p^}} , (6) 

for B s -> D*+D*-. In above equations, g = ^V cb V* s , fj* ~ 200MeV, f±,V ~ 0.6 and 
A\,Ai — 0.25, according to Ref. ||. The matrix element for the decay B s — ► 7(51^)7(921/) 
are parametrized as, 
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A = ag (RiS^ + iR 2 P^) , 

Pfiv = e (Uva/ gg"g 2 j (7) 

where i?i and _R 2 are the yet-to-be-determined hadronic matrix elements, and the partial 
width is then given by, 

r(5 s -77) = ^m|(|i? 1 | 2 + |i? 2 | 2 ) . (8) 
Since T(B S ) ~5x 10- 10 Me\/, we have, for \V cb \ = 0.04, 

1-7 f\ ^1 1 2 i i -^2 



In re-calculating the processes, for the DfD~ intermediate state we confirm the result 
of Ref. |U, Eq. (|l|). For the D* + D*~ intermediate state, we have, 



lmR 1 = -^^{[(l - 4m) + y (4m) - <fo]\ 
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+m[(l - 12m + 48m 2 ) + (f> v {-2 - 8m) + (p 2 v (l - 4m)] My^)} 
^^{[(1 - 5m) + y (2m) + $(-1 + 3m)]A 22 



8m 



+m[(l - 10m + 32m 2 ) + y (-2 + 4m) + <$(1 - 2m)] M^^)} , 

ImR 2 = ^-{[(-1 + 12m) + y (2 - 8m) + <^>(-l - 4m)]A 22 
lorn 

+ [(4m - 32m 2 ) + y (8m - 32m 2 ) + <^(4m)] log(^4)} . (10) 

The calculation made in ref. || corresponding to taking <py = in above equations. We 
find a disagreement on the first term inside the square bracket in the coefficient of A 2 . As 
a consequence, the CP conserving branching ratio becomes much smaller, BTZ(B S — > 2 / ~f) ~ 
0.85 x 10~ 8 . Meanwhile we reconfirm the numerical result on the CP violating branching 
ratio of ref. M. The maximal value of B1Z(B S — > 27) comes from the case <py — 1. In this 
situation, our calculation is equivalent to taking the W + W~ , y vertics (propagators are in 
the unitary gauge). Hence we list some intermediate steps in the calculation here since it 
may be useful elsewhere. The s partial-wave projection of the 2 — > 2 scattering amplitudes 
are found to be, 

PT.(d+d; - 27) = g(i - P 2 ) ln(^|) , (11) 

and, 

pT++(D?Dr - 2 7 ) = g(l + P) 2 In (l±E\ , 



/5T s -(^r-2 7 ) = |(l-/5) 2 ln(i±|j , 

CTW - 27) = - /? 2 ) 2 In (I±|) . (12) 

In above the gauge invariant factor ei • e2&i ■ ki — e\ ■ k 2 e 2 • k\ (= s/2) has been taken away, and 
the photon helicity is taken as (+,+)• Once again we find that the D* + D*~ intermediate state 
contribution to the decay branching ratio of B s — > 77 quite small, though it is comparable 
in magnitude to the short distance contributions, 

BK CP -even{B s -> 2D* S -> 2 7 ) ~ 1.49 x lO" 7 , Bn C P-odd(B s -> 2D S * -> 2 7 ) ~ 0.44 x 10" 7 . 

(13) 

From the above results we find that the decay branching ratio of B s — > 77 will not exceed 
a few times 10~ 7 provide that the real part contribution does not change the above estimtate 
in order of magnitude. It is natural to expect that the real part contribution is comparable 
in order of magnitude to the absorptive part. Therefore, if the future experiments reveal a 
large branching ratio of B s — > 2 7 upto a few times 10~ 6 it is very likely that it is a signal of 
new physics. 

Now we discuss the real part contribution of the process depicted in fig. [I]. The real part 
contribution can be calculated using dispersion relation. When performing the dispersion 
integral it is necessary to extract the e\ ■ e 2 k\ ■ k 2 — €\- k 2 e 2 ■ k\ part out of the dispersive 
integral |10 since it is the remaining part to be considered as the form-factor. [] We get, 



A = I I"" ' mA ds> . (14) 

TT J 4m 2 S' — S — ie 

The D s loop in fig. 1 contains no ultra-violet divergence and it is straightforward to 
evaluate Eq. (|T4"D, we find, 

BU{B S -> 2D S -> 77) ~ 4.2 x 10" 8 . (15) 

For the D* contribution the diagram fig. 1 is ultra-violet divergent and needs to be renor- 
malized in the standard perturbation calculation. In the dispersive approach the divergence 
manifests itself by the fact that the dispersive integral in Eq. ( JT4] ) needs one subtraction. 
The divergence can also be rescued by introducing form-factors in the Feynman vertex to 
regulate the high energy behaviour of the Feynman diagram. However, it is equivalent in 
principle to regulate Eq. ( |I4] ) by using an un-subtracted dispersion relation with truncated 
integrand, 

1 r A2 ImA 

A = - / ds' . (16 

TT J 4m 2 S' — S — ie 

Of course, the ultra-violet divergence emerged from the effective Lagrangian approach soley 
indicates that the effective Lagrangian approach gives a bad high energy behaviour to the 



1 For earlier referrences on the usage of dispersion relation under such a situation, see ref. |1C,11| 



5 



rescattering amplitude. For example, the T s ++,++ amplitude of Eq. (|12[) , behaves as ~ - ln(s) 
as s — > oo and does not satisfy Eq. fll4]) . However, from Regge theory we know that, at high 
energies, the s-wave scattering amplitudes behave as T s (s) ~ s 00 *' 1 and T s (s) ~ s" 1 ^ -1 , for 
the D s and .D* exchanges, respectively. The intercept parameters of the Regge trajectory are 
ao s — —1.5 and ao s * — —1- It is clear that the dispersion relation, Eq. ([HP, is convergent. 
The rapid fall-off behaviour of the Regge amplitude with respect to s implies that the 
magnitude of the real part contribution obtained in the effective Lagrangian approach and 
via Eq. (^) be well overestimated. In fact, Regge behaviour should already dominate at 
s = M\ s [JEJ. So the result obtained from Eq. ( |To| ) may at best be interpreted as an upper 
bound. Since the high energy contribution can not be important to the dispersive integral 
in Eq. ([14]), we take A ~ 2M^ s as an educative estimation. We find that the inclusion of the 
real part contribution does not alter the order of magnitude of the absorptive contribution. 
For example, when taking 0y = 1, we have, 

B1l(B s -> 2D* S -> 2 7 ) ~ 3.7 x 1(T 7 , A = 2M 2 Bs , (17) 

to be compared with Eq. ([13]). 

To conclude, we have re-evaluated the B s — ► 27 decay process through D+D~ and 
D* + D*~ intermediate states' rescatterings. We have estimated the absorptive part contri- 
bution to the decay branching ratio in the effective Lagrangian approach, and also estimated 
the real part contribution through a cut-off regulated dispersion relation. We argue that the 
decay branching ratio will not exceed a few times 10~ 7 . Therefore if the furture experiments 
reveal a large branching ratio of a few times 10~ 6 it must come from new physics effects. 
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